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We study the fluid-like dynamics of eigenvalues of the Wilson operator in 
the context of the order-disorder (Durhuus-Olesen) transition in large N c Yang- 
Mills theory. We link the universal behavior at the closure of the gap found by 
Narayanan and Neuberger to the phenomenon of spectral shock waves in the com- 
plex Burgers equation, where the role of viscosity is played by 1/N C . Next, we 
explain the relation between the universal behavior of eigenvalues and certain 
class of random matrix models. Finally, we conlude the discussion of universal- 
ity by recalling exact analogies between Yang-Mills theories at large N c and the 
so-called diffraction catastrophes. 

PACS numbers: PACS numbers come here 

1. Introduction 

Many efforts continue to be devoted to the study of QCD in the limit of a large 
number of colors, after the initial suggestion by t'Hooft (T). This, in part, is due to 
the general belief that the large N c limit captures the essence of confinement, one 
of the most elusive of QCD properties. At the same time the theory simplifies con- 
siderably in the large N c limit: fluctuations die out and the measure of integration 
over field configurations in the partition function becomes localized at one par- 
ticular configuration, making the large N c limit akin to a classical approximation 
0. 
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Many results have been obtained in the simple case of 2 dimensions. Then 
Yang-Mills theory translates into a large N c matrix model, where the size of the 
unitary matrix is identified with the number of colors. More specifically, the basic 
observable that one considers is the Wilson loop along a (simple) curve C 

W[A] = Pe'&'V^ (i) 

where A M = A^T a , with T a the generators of S\J(N C ) in some representation. In 
the fundamental representation, W[A] is an N c x N c unitary matrix with unit de- 
terminant. Its eigenvalues are of the form A = exp(/#) and can be associated with 
points on the unit circle. After averaging over the gauge field configurations, with 
the usual Yang-Mills measure, one finds that W =< W[A] > depends in fact only 
on the area enclosed by C, to within a normalization 0. It is convenient to 
measure the area in units of the t'Hooft coupling g 2 N c , i.e., we set r ~ g 2 N c J{. In 
the limit N c — > oo, the eigenvalues are distributed on the unit circle according to 
an average density p(6, t). 

The typical behavior of the distribution of eigenvalues as a function of the 
area is displayed in Fig. [T] One observes that for small loops (which probe short 
distance, perturbative phyics), the spectrum does not cover the whole unit circle, 
but exhibits a gap; in contrast, for very large loops (which probe long distance, 
nonperturbative physics) the spectrum covers uniformly the unit circle (gapless 
phase). This behaviour of the spectrum agrees with the order (gapped)-disorder 
(gapless) transition, proposed by Durhuus and Olesen [3 ] and based on the explicit 
solution of the corresponding Makeenko-Migdal equations in 2 dimensions. 
Surprisingly, a similar critical behavior has been observed also in d = 3 dimensions 
and conjectured to hold in d = 4 large N c Yang-Mills theory [5], suggesting a 
universal behavior (see the lectures by Neuberger and Narayanan in these volume, 
to which we also refer for a discussion of the subtleties of the regularization of the 
Wilson loops). 

This universality conjectured by Narayanan and Neuberger is comforted by 
simple schematic matrix models, in particular that proposed by Janik and Wiec- 
zorek [6], hereafter JW model. The model stems from the general construction of 
multiplicative free evolution Q, where increments are mutually free in the sense 
of Voiculescu OEI. The unitary realization in the JW model corresponds to ma- 
trix value unitary random walk, where the evolution operator is the ordered string 
of consecutive multiplications of infinitely large unitary matrices 




where Uk = exp/l yt/KHk), with a hermitian random matrix, drawn from a 
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Fig. 1 . The evolution of the spectral density as a function of the area J{ of the average 
Wilson loop (d=2). The plots of p(9,^i) at the bottom are taken from the simulations of 
Ref. 0. 



Gaussian probability distribution P(H) of the form 

P(H) ~ e - OTry(ff) <TfTr#> - <^Tr// 2 > = m 2 . (3) 

The model is a random matrix generalization of the multiplicative random walk 
performed in K steps during "time" t. In the continuum limit K —> oo, the model 
is exactly solvable. The solution for the spectral density coincides exactly with 
that of the two-dimensional QCD, provided one identifies t with the area of the 
Wilson loop, modulo a normalization ifTUl . This model offers a neat picture for 
the multiplicative evolution: at t = 0, the spectrum of W is localized at A = 1. As t 
increases, the spectrum starts to spread symmetrically along the unit circle towards 
the point A - - 1 , reaching this point and closing the gap at finite time. Further evo- 
lution corresponds to further spreading of eigenvalues around the circle, resulting 
finally in a uniform distribution (see Fig.[T]). 

Neuberger and Narayanan Q have observed, that large Af c Yang-Mills lattice 
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simulations in d = 2 and d = 3 demonstrate the same critical scaling at the closure 
of the gap as in the JW model and have conjectured that this model establishes a 
universality class for d = 4 large N c Yang-Mills theory as well. In their simula- 
tions, Narayanan and Neuberger [5] did not calculate the spectral density directly, 
but rather the average characteristic polynomial 

Q Nc (zJ) = {det(z-W(t)))- (4) 

As we shall see later, this contains the same information as the spectral density 
when N c — > oo. Narayanan and Neuberger performed simulations at finite N c 
and obtained evidence that the crossover region between the gapped and gapless 
regimes is becoming infinitely thin as N c — > oo. 



2. Spectral density and resolvent 

The object at the heart of our discussion will be the average density of eigen- 
values p(9, t), defined so that the number of eigenvalues of the Wilson operator in 
the interval [9, 9 + dd], after averaging over the gauge field configurations loops of 
a given area J{ ~ r, is p(9, r)d9. 

2.1. The spectral density and its moments 
The spectral density p{9, r) is not available in analytic form, but its moments 

X+7T 
d9e in0 p(9,T) (5) 
n 

are. An explicit, compact form for these moments is given in Ref . |[T4ll in terms of 
an integral representation 

w n (r) = + l/zf cxp(-nT(z + 1/2)) 

n J 2m 

= ^L ( 1 „_ 1) (?iT)exp(-?iT/2) (6) 

where the representation of Laguerre polynomials, used in the second line, allows 
connection to results known already 25 years ago |f3j |TT| . 

The Durhuus-Olesen transition can be seen by studying the asymptotic behav- 
ior of these Laguerre polynomials, using a saddle point analysis of their integral 
representation lfT2l[T3"l . The result is surprising: for a loop area below the critical 
value t c = 4, the moments oscillate and decay like n~ 3/2 , while for r > r c the mo- 
ments decay exponentially with n, modulo similar power behavior. Both regimes 
are separated by a double scaling limit. 
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Let us quote here the values of the first couple of moments. The normalization 
of the spectral density is given by 



w = f dGp{6) = \. (7) 

The first moment expresses the area law obeyed by the average of the Wilson 
operator 

m = e~ r/2 . (8) 

More generally, since p is real, w* = w_„, and since the moments, as given by the 
formula above are real, we have (n > 0) w_„ = w n . Thus we can express p as 
follows 

p(0) = y Z WmQ ~ im9 



2n 

m=-oo 



1 ( + °° 

= — 1 + / 2w„ cos(«^) 

V n=\ t 



(9) 



Note that, for r = 0, w n = 1 for all «, and 



1 +oo 

m=-oo 



= J] 5(6l + 2«7r) 



n=-oo 

= 5(0) (10) 
where we used Poisson summation formula, and in the last line the restriction 

-7T <6 <7T. 

2.2. The resolvent 
To the spectral density one may associate a resolvent, defined as 

G(z) = £d9f^. (11) 

For \z\ > 1 we can expand the integrand in a power series of e l6 /z and get 

CO 

zG(z) -1 = 2^ = /(z) ' (12) 
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where the function / is the same as that introduced in Ref. [6]. By setting z = e' a , 
we can also write 

oo 

/(^e'Vfr.e-'" (13) 

n=\ 

with Im a < to guarantee the convergence of the sum. 

Now we introduce a new function that will play a central role in our discussion. 
We set |H1: 

F(z = e ia ) = ilf(z) + l) = i(zG(z) - I 



= ^ + |>»e-^j, (14) 

where in the last line we assume Ima < 0. The imaginary part of F yields the 
spectral density (8 real) 

p(0) = -lmF(z = e w ). (15) 
We can also write (a and 9 real) 

n JW) {e^T^-2/' (16) 



so that 

1 C +n I a — G\ 

F(z = z ia ) = - J d6p{0) cot y— J . (17) 

At this point, we shall make an abuse of notation and write F(a) for F(z = e' a ), 
and furthermore we shall allow a to be complex. One can then write, for real a 

-F(a) - Hp(a) + ip(a), (18) 

7T 



where Hp(a) is the Hilbert transform 

1 r +n (a-9\ 
Hp(a) = — P.V.J dSptf) «*(—). (19) 

Note that there is a choice of sign for the imaginary part, related on where we 
decide to view F as analytic (lower or upper plane), that is, on how we take the 
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limit from complex a to real a. The present choice corresponds to choosing F to 
be analytic in the lower half-plane. 

The following properties of the Hilbert transform will be useful 



He ia = i-P.V. J* *d0e*cotp^) 



= -ie ia . (20) 

ztt j_ n \ z I 

It follows that 

// cos a = sin or H sin a = - cos a (21) 

We have also 

H(H(f)) = -f 

H{fg) = fH(g) + gH(f) + H[H{f)H{g)], (22) 

from which it follows in particular that 

H[H(p)p] = ~((H(p)f - p 2 ) . (23) 



3. Complex Burgers equation and its solution with characteristics 

The usefulness of the function F that we have introduced comes from the fact 
that it satisfies a simple equation, the complex Burgers equation |[T4l [3l [TTl : 

d T F + Fd F - 0. (24) 

This equation is analogous to the real Burgers equation of fluid dynamics lfT31 
(with r playing the role of time, 9 that of a coordinate, and Fq of a velocity field). 
The complex Burgers equation is omnipresent in Free Random Variables calcu- 
lus flU. It also appears frequently as one dimensional models for quasi-geostrophic 
equations, describing e.g. the dynamics of the mixture of cold and hot air and the 
fronts between them ifToll . Here, we shall take advantage of the abundant mathe- 
matical studies of the complex Burgers equation to analyze the flow of eigenvalues 
of Wilson loop operators. We shall in particular use the crucial fact that the com- 
plex Burgers equation may allow for the formation of shocks. 

A simple derivation of the complex Burgers equation will be given in Sect. [5] 
below. In this section we shall analyze general properties of its solutions, using the 
method of complex characteristics. 
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Fig. 2. Characteristic lines for t < 4, a — £ + \ cot |, with a - 8 real on the vertical 
axis, and t on the horizontal axis. The various straight lines are given by the characteristic 
equation ( 33 i in which £ = x is treated as a parameter. The horizontal characteristic 



corresponds to x — n, the first one, almost vertical, to x = 0.1. 



3.1. Characteristics 

The Burgers equation 

d T F + Fd a F = 0. (25) 
admits the following solution in terms of characteristics 

F(r,a) =F (f(T,a)), F (a) = F(t = 0, a), a = Z + TF (&. (26) 
The initial condition Fo corresponding to a spectral density peaked at 9 = 0, 

Po(A) - <5(0), (27) 

is 

a 



1 

fo(ff) = 2 cot 2 ■ 
The characteristics are therefore given by 

ar = f+Icot|, 



(28) 



(29) 
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Fig. 3. The characteristic function 6{x,y) (left) and T](x,y) (right) for t = 0.5, cut by the 
plane — 1.38 and rj — 0, respectively. 



with a complex. Once g(a, t) is known, F(a, t) can be obtained as F(a, t) = 
(«-^)/t. In particular, for real a, we have ImF = -(l/r)Im£, so that p = -y/{rn), 
with y = Im£. Alternatively, one may look for F as the solution of the implicit 
equation 

F = icot(£ + TF). (30) 

We shall set £ = * + (y, and a = + irj. Then, by taking the real and imaginary 
parts of the characteristic equation, we get 

sin* t sin* 

9 - x + re- y — - x + 



1 - 2e y cos x + e 2 >' 2 coshy - cos * 

t 1 - e 2v t sinhy 

r i = y + ^-y — ^r~v ^ = y~^, — u ■ ( 31 ) 

2 1 - 2e> cos x + & l y 2 cosh y - cos x 

The characteristics form a family of straight lines in the (9, if) plane, 9{t), t](t), 
that depend on the paramaters x and y (the values of 9 and rj at t = 0). Since the 
functions 0(jc, y) and rj(x, y) are the real and imaginary parts of an anayltic function, 



a{%) given by Eq. (29), they satisfy the Cauchy Riemann conditions: 



89 A cos*coshy-l drj 89 sinxsinhj drj 

8x 2 (cosh y - cos*) 2 dy' 8y 2(coshy - cos*) 2 dx 

(32) 

An example of characteristic functions 9(x,y) and rj(x,y) is displayed in Fig. [3] 
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As a simpler illustration, we note that the characteristics that start at y = 0, and 
< x < 7r, remain in the plane 77 = as t varies. They are plotted in Fig. [2| As 
easily seen from eqs. pTj ), when y = 0, rj = 0, and 



T X 

= x+ - cot-. (33) 
2 2 



The envelope of this family of lines, also drawn in Fig. [2] is given by 9 C = x c + 



(t/2) cot(;e c /2), where ;e c is obtained from the equation (see eq. (43 1 below) 



d x e = o = l = . (34) 

4sin 2 (x/2) 

Let us emphasize that these characteristics in the plane 77 = are not enough to 
construct the solution of the complex Burgers equation. In fact, by restricting the 
starting coordinate to be real (y = 0), we have produced a set of characteristics 
that cross each other. This set of characteristics is "unstable": As soon as a small 
imaginary part is present (i.e., y + 0), the characteristics move away from the plane 
77 = 0. 



3.2. Graphical solution 

To obtain the solution of the Burgers equation from the characteristics one 
needs to identify the characteristic (i.e. determine its origin £ in the plane r = 0) 
that goes through the point (9, 77) at time t. In other word, one needs inverting 
the relation a(£, r) — > %{a, t), and inserting this value of £ in Fq(£). This can 
be done graphically by first drawing the lines of constant 6 and constant 77. We 
are interested mostly in the values of F for real angles a, that is for 77 = 0_. 
Fig. [4]provides an example for t = 0.5. This figure could be used as a basis for 
a "graphical" solution of the Burgers equation: each intersection point in Fig. [4] 
gives the coordinates x, y of the origin in the plane t = of a characteristic going 
through the point of coordinates 6, 77 = (where 9 can be read on the corresponding 
level line). We shall examine later the cases where this construction fails. 



3.3. The curves of constant rj - 

Because p ~ y, it is instructive to analyze the shape of the curves 77 = in the 
(x,y) plane. Indeed this shape will be similar to that of the function p(9), with 9 
related to x as described in the caption of Fig.|4] Some curves 77 = are displayed 
in Fig. [5] These are determined by the implicit equation 
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Fig. 4. The set of curves represent level lines of constant 9 in the x,y plane. These lines 
are intersected by the line of constant 77 = (77 < "inside", and r\ > "outside"). 
Each intersection point represents the origin (x,y) of a characteristic arriving at point 
(9,rj = 0) in time t. Here t = 0.5. The various lines corresponds, from left to right, 
to 6 = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8, 1, 1.2, 1.3, 1.35, 1.38, 1.384, 1.3845, and the line 
crossing orthogonally the lines of constant 9 is the line 77 = for the same value of t. Note 
that there is a solution (i.e. an intersection point) only for x ~ 0.7. This point is the sin- 
gular point associated to the edge of the spectrum, given by £ e = x c — 2Arcsin( VO.5/2) 
0.722734, corresponding to 8 C = 0.7227234+ V0.5(l -0.5/4) ^ 1.38417. Note that as one 
approaches this singular point, the landscape becomes completely fiat, i.e., d x 6 = d y 9 = 0. 
Note also that the curve of constant r\ — is essentially the spectral density p(6, t), with 
9 determined form the intersection points (this follows from the fact that p = —y/(nr) 
for negative y, and the fact that the set of lines in the figure is the mirror image of the 
corresponding set fory < 0). 



It is easy to show that the curves rj = intersect the y-axis with a vanishing 
slope. Let us focus on their behavior near the x axis. For small y we have 

t y t 

y = -- , cosx c = l--. (36) 

2 1 - cos x 2 



Thus, as long as r < 4, there is a value x c of x at which the curve rj - intersects 
the ^-axis. As we shall verify later, this point is the singular point associated with 
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Fig. 5. The lines of constant 77 = in the x,y plane. The different curves correspond 
to t — 0.5, 1,2,4,4.1,5,7. These lines yield also the density p(r,x). In order to get the 
density p(r, 6) we need the mapping between x and 0, which can be read on a figure like 
Fig.0 



the edge of the spectrum. In the vicinity of this point, the curve 77 = does not 
depend on y, that is, it has infinite slope in the (x,y) plane. 

For r > 4, the curve intersects the axis x = n at a point y s solution of 

y s = ^tanh|. (37) 

y s is a growing function of t. For large t, y s - ±t/2, and one can easily construct 
the solution. Indeed we have then £ = n + iy s ~ n ± it/2, so that 

F O (0 = \ cot n _l^Jl = ± L tanh I (38) 



and the density is 



pfr,T) = — tanh^. (39) 

27T 2 



For r — > 4 + 0, the equation for y s becomes 



Tiy y 3 \ 2 
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so that 

P(7T,T)- i-V3(T-4). (41) 

Thus the derivative dp/dr is singular at r = 4. 

3.4. 77ze caustics 

The construction of the solution of the Burgers equation from the character- 
istics is possible as long as the mapping between a and £ is one-to-one, that is, 
as long as da/dg + 0. When da/dg = a singularity develops. The equation 
which determines the location of the singularities is also that which determines the 
envelope of the characteristics, the so-called caustics of optics. Let £ c (r) be the 
location of the singularity. We have 

F'otfc) = - . I, - =~. (42) 
4sin 2 (&/2) r 

The equation for the caustics is then given by 

a(r) = &(t) + T - cot (43) 

By setting £ c = x c + ry c , we transform the singularity condition into two equations 
for x c and y c : 

sinh(y c /2) cos(x c /2) = 

cosh(j c /2) sin(x c /2) = ± Vr/2. (44) 

These equations are equivalent to the equations d x 9 = d y 9 = (we used the 
Cauchy-Riemann conditions). The first equation implies that y c = unless x c = 
n. Consider the first possibility, i.e., y c - 0. The second equation then yields 
sin(x c /2) = ± ^rll which is possible only if t < 4. One concludes therefore that 
if t > 4, y c + and x c = n. Consider now the second possibility, x c = n. In this 
case, the second equation yields cosh(y c /2) = ± which is possible only if 

t > 4. Therefore if r < 4, x c + n and y c - 0. In summary, for r < 4, the caustic 
lies in the plane r\ - 0, while for t > 4 it lies in the plane 9 - n. 
For r < 4, the value of x c is given by 

x c - 2 arcsin( yfr/2) - arccos(l - t/2). (45) 

The equation of the caustics in the 9, rj plane, is given by 

9 C = 2 arcsin( Vr/2) + Vt(1 - t/4), ?7 c - 0. (46) 
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This is the curve plotted in Fig[2| The value 9 C corresponds also to the edge of the 
spectrum when t < 4. When the gap closes, t — > 4, 9 C — > £ c = nr. 
For r > 4, j c is given by 

j c = 2argcosh a/t/2, (47) 

and (for y > 0) 

T] c = 2 argcosh Vr/2 - Vt(t/4 - 1), (48) 

with -rj c also solution (for y < 0), describing a symmetric branch of the caustic. 
Note that as r — > oo, ?7 C . -t/2. 



3.5. Solution of Burgers equation in the vicinity of the caustics 

In the vicinity of the caustics one can construct the solution of the Burgers 
equation analytically. This is because one can then easily invert the relation be- 
tween £ and a. One has, quite generally, 

F (g) = F (&) + (f - &)*o(&) + \(g ~ ZcfF'Jitc) + i(f - ZcfFo'igc) + ■■■, 



so that 



(49) 



a = a c + U% - tcfF'^c) + X -M ~ gcfF'o'tfc) + • • ■ (50) 
2 o 

For t < 4, we have £ c = 2 arcsin( y/r/2) and 

tF (&) - V(1-t/4), 



For t < 4, one can ignore the cubic term. One gets then 

a = a c + (g - £ c ) 2 VI/t- 1/4, (52) 
which is easily inverted to yield 

f - & = — T77 V» - ffc- (53) 

s s (1/t- 1/4) 1 / 4 v 
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We have then 



F(a, r) = F (f (or, r) - F (£ c ) + (g - &)F' (&) 



V t(1/t- 1/4) 1 / 4 V 



The spectral density can be deduced from the imaginary part of F, or equivalently 
,Fo-We have 

p(a) = -LnF(o - i0+). (55) 

It follows that the spectral density vanishes for a > a c (confirming the interpreta- 
tion of a c as the edge of the spectrum). For a < a c , 

P(«) - n/ 1 , ; , u;4 y/ac-a . (56) 
/rr(l/T - 1/4) 1 ' 4 

When t = 4, the second derivative vanishes, a c = n = % c , and we have 

a = n-^-(t-% c ) 3 , a-n = ~(t-n)\ (57) 
3t 3t 

so that 

f - 7T a (3T) 1/3 e i7r/3 (o: - tt) 1/3 . (58) 
It follows that the spectral density is given by (with a real and a < ri) 



m 



1/3 

p^.^sin^-ay/^i-^j (n-af'\ (59) 



As r grows beyond r = 4, the real caustic splits into two complex ones, mov- 
ing in opposite directions along the 77 axis. At this point we shall not pursue our 
analysis with complex characteristics, but go through a specific study at 9 = n 
using a different approach. 

3.6. Specific study at 9 — n 

In this section, we study the solution for 9 = n, using an approach which will 
allow us to make contact with the work of Neuberger fP71IT8l . 



16 



BN'Zakopane'09a printed on November 19, 2009 



3.6.1. Solution with a real Burgers equation 

Starting from the complex Burgers equation, for complex a, 

d T F + l -d a F 2 = 0, (60) 

we restrict a to a = n + irj - n - iy (where y is the variable introduced by Neu- 
berger). Choosing y > 0, we get r\ < so that we are in the right domain of 
analyticity of F. In fact, F is analytic except possibly at y = (where it develops a 
discontinuity when the gap closes; this discontinuity is proportional to the spectral 
density at a = n). To avoid confusion (in all this discussion x and y will not have 
the same meaning as before in these lectures) we use Neuberger's notation and set 

0(y) = iF{a = tt- iy), drfry, r) + -d y <f> 2 (y, r) - 0, (61) 
with the boundary condition 

0(y,T = O) = 0oOO = ~tanh;|. (62) 
Note that since F{n + ie) = ±in(p(n), we have 

p(/r, t) = — [0(O_, r) - 0(0+ , t)] . (63) 



3.6.2. Characteristics of the real Burgers equation 

The characteristics are given by 

T X 

y = x + T<f> (x) = x- 2 tanh 2' ( 64 ) 

and the solution is given in terms of them by 

<f>(y, t) = <po(x(y, t)), x(y, t) = y - T<p (x). (65) 

The velocity field <f>(y, t) will eventually become infinitely steep at some points y c . 
To determine these points we calculate the derivative d y (p(y, r): 

d(f> d(f>o dx d(po 1 
<9y dy dx dy/dx' 

where we use the fact that y and x are related by the characteristic equation. It 
follows from this equation that 
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Fig. 6. The function f T (x) plotted as a function of x for the values r = 3,4,5. 



Since d<f>o/dx in Eq. (66 1 is finite for all x, an infinite slope in the velocity field 
will develop when dy/dx = 0, that is, for values of x c such that 



1 +T 



d<f>o 
dx 



0. 



(68) 



It is convenient to set 



Mx) = <po(x) + 



(69) 



The characteristic equation ( 64 1 reads then y = rf T (x) = 0, and the location of the 
singularities is given by f'(x c ) = 0. The relation between x and y can be inverted 
if the solution of the equation f T (x) = y/r is unique. The plot in Fig. [6] indicates 
that this happens when r < 4. Indeed the derivative of f T (x) is given by 



i 



i 



r 4cosh 2 (x/2) 



/r'(0) 



1 1 



(70) 



so that the derivative in x = is positive as long as r < 4. In this regime, there is a 
unique solution x for each y. When t = 4, /'(0) = 0, i.e., x = = x c and the slope 
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at y = of <p(y, 4) is infinite. This is the preshock corresponding physically to the 
closure of the gap and the beginning of the build up of the spectral density at 6 = n. 
For t > 4 the solution becomes multi-valued (takes an S-shape). The function 
cp becomes then discontinuous, its discontinuity giving the density, according to 



Eq. (63 1 



3.6.3. Solution of the real Burgers for large r 

For large r, the solutions of the equation f T (x) = are approximately given by 

T T T T 

;q = -tanh- (x > 0), ;q = --tanh- (x < 0). (71) 

The solution of the characteristic equation y = t/ t (x) are shifted linearly with 
respect to these values, when y is small, that is, ±xi — > ±x\ + y. It follows that 
<P(y> t) - (y — xi)/r becomes independent of y in the vicinity of y = 0, except for a 
jump that depends on the sign of y. We have 

<f>(y < 0,t) * i tanh T -, </>(y > 0, r) - -i tanh T -. (72) 



It follows that 



1 r 

P(tt,t) - —tanh-, (73) 

Z7T 2 



in agreement with Eq. (39 1. 



3.6.4. Introducing a small viscosity 

Anticipating on the discussion in the next sections, it is interesting to consider 
the viscid Burgers equation for <f>(y, r) 

d T cp + l -d y <p 2 = vd 2 y <(>, (74) 

where v plays the role of a viscosity (we shall see later that v = 1 /2N C ). This 
equation can be solved with the so-called Cole-Hopf transform 

<f>(y,T) = -2vd y In K(y,r), (75) 

where K satisfies the diffusion equation 

dK _ v d 2 K 
8t dy 2 

Let Kq(x, t) be the solution that reduces to 6{x) at r = 0: 

Kq(x,t) = — =e-^. (77) 
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Then the solution of the Burgers equation reads 



We have 



K(y.rt=— — dxe-^e^i>(«y«. (78) 
yArtvr J-oo 



I du(f>Q(u) = -In cosh -, (79) 
Jo 2 



and we can write 



1 



X+oo 
CO 



1 



K(y,r) = — ^e~fc | dxe~tv v (-x\ (80) 



with 



2t t 2 192 \t 4/ 2 t v ' 

The expansion on the r.h.s. allows us to recover the Pearcey integral and study the 
vicinity of the critical point (near t = 4) and the scaling with (small) viscosity. We 
have indeed 

X+oo /^\^ 4 r-+°° 

dxe~T v V(x) = 4\±\ v l/4| dMe -« 4 -a« 2 + f« (g2) 



x = 4(^l v 1/4 M , (83) 



with 

s 2 

and the scaling variables related to r and y are, respectively 



4. The gapless phase and the inverse spectral cascade 

In this section, we shall illustrate a particular feature of the disordered (gap- 
less) phase. Consider the large r uniform solution, and a small perturbation of the 
spectral density of the form 

p(0,r o ) = t"-(1 +26 cos 0) (85) 
2n 
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with to » 1 . Note that this form of the density results from truncating the general 
expansion (|9]) at the first moment, and set w\ - e. We wish to solve the Burgers 
equation with p{9, To) as the initial condition. The characteristics are given by 

= f + (T-T O )Fo(f), (86) 



and the function Fo(£) corresponding to the initial condition can be read off Eq. ( 14 1: 

Fo<& = \ (l + 2ee"*) . (87) 
A singularity occurs for f = g c , with g c solution of 

e& = -e(r - T ). (88) 

At g c , we have 

Fo(^ c ) = \ - - J —, = — , = (89) 

2 T-T() T-T T-T 

We may now proceed and determine the solution in the vicinity of the singu- 
larity, as we did in Sect. [3] In the vicinity of the singularity we have 

6 = e c + l -{t-z c f. (90) 



The equation for the singularity, Eq. ( 88 1, has two solutions, depending on whether 
t > to or t < To. Let us consider these solutions in turn. 

If T > T(), 

t; c - n - i In £{t - To) + 2nn, (91) 

and 

e e = n-i{l- + In £(t - T )) (92) 

When t is near To, T > to, the singularity is at 9 C ~ n + ioo. As t keeps increasing, 
6 C remains complex unless e is too large (i.e., unless e > 1/2). So, for small e, as 
one propagates forward in time, B c remains complex and return to +ioo as t — > oo. 
Note that at large t, Fq(^ c ) ~ i/2 and the density p ~ 1 /2k. 
The situation is different for x < to. Then 

& = -i In e(T - t) + 2nn (93) 

and 

Q c = -i 1 1 + ^(T - t) + In e(T - t)J (94) 
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Again the singularity is at +ioo when t = To, but as t decreases, it moves towards 
the real axis and reaches it in a finite time t* given by 

0= l + i(T -T*) + ln6(T0-T*) (95) 

At this point we have the usual blow up, and the solution ceases to exist. 

This simple calculation shows that while the complex Burgers equation pro- 
vides a unique solution that connects the singular density at t = (p(9) = 6(d)) and 
the uniform density at r — > oo (p = l/2n), this solution is "unstable" to backward 
motion for initial conditions that deviate only slightly form the generic solution 
of the Burgers equation. It is tempting to speculate that this particular feature is 
related to turbulent aspects of the disordered phase |[T9l . and indeed the Burgers 
equation has often been used as a toy model to study such turbulent behavior. As 
a further remark, note that as one propagates forward in time, the higher moments 
of the spectral density are damped, leaving eventually only the uniform density at 
large time. This process is reminiscent of the inverse cascade, a feature of two- 
dimensional turbulence that can also be studied with the Burgers equation (for a 
pedagogical illustration of the phenomenon, see BUlO . 

5. Dyson fluid 

In this section we would like to demystify somewhat the omnipresence of in- 
viscid and viscid Burgers equations in the analysis of the closure of the spectral gap 
of Wilson loops. We will see that these equations originate from either algebraic 
or geometric random walks, where however the standard independent scalar incre- 
ment is now replaced by its matrix-valued analogue. Actually, the idea of matrix- 
valued random walks appears already in pioneering papers on random matrix the- 
ories. In a classical paper EH . Dyson showed that the distribution of eigenvalues 
of a random matrix could be interpreted as the result of a random walk performed 
independently by each of the matrix elements. The equilibrium distribution yields 
the so-called "Coulomb gas" picture, with the eigenvalues identified to charged 
point particles repelling each other according to two-dimensional Coulomb law. 
For matrices of large sizes, this correctly describes the bulk properties of the spec- 
trum [22 J. In his original work, Dyson introduced a restoring force preventing 
the eigenvalues to spread for ever as time goes. This is what allowed him to find 
an equilibrium solution corresponding to the random ensemble considered with a 
chosen variance (related to the restoring force). We do not need to take this force 
into account here, since the phenomenon we are after is a non-equilibrium phe- 
nomenon. In particular, in the picture that we shall develop, and that we may refer 
to as "Dyson fluid", the edge of the spectrum appears as the precursor of a shock 
wave, and its universal properties follow from a simple analysis of the Burgers 
equation that was developed in other contexts E3ll . This observation allows us to 
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link familiar results of random matrix theory to universal properties of the solution 
of the Burgers equation in the vicinity of a shock EM . 

5.1. Additive random walk of large matrices 

We start from the simpler case of additive random walk of N x N hermitian 
matrices H with complex entries. The explicit realization of the random walk is 
provided by the following construction: In the time step 8t, Hij — > Hjj + 5Hjj, with 
(SHij) = 0, and ((6Hij) 2 ) = (1 + 8ij)6t, so at each time step, the increment of the 
matrix elements follows from a Gaussian distribution with a variance proportional 
to 5t. The initial condition on the random walk is that at time t - 0, all the 
matrix elements vanish. The crucial difference between the scalar and matricial 
random walk is visible when we switch from matrix elements to eigenvalues of 
H, which we denote by a,. The random walk of the eigenvalues has the following 
characteristics EH 

{5x i ) = E{x i )5t, {{6xif) = 5t, (96) 
where the "Coulomb force" 

^-Zijfh) < 97 > 

originates in the Jacobian A of the transformation from the matrix elements to the 
eigenvalues, A = Yli<j( x i ~ x j) 2 - Note that now random walkers interact with 
each other, due to the "electric field" E(xj - Xj). This interaction may a priori in- 
troduce non-linearities in the corresponding Smoluchowski-Fokker-Planck (SFP) 
equations. Indeed, this is the case, as we demonstrate below. 

Using standard arguments, we see that the joint probability P{x\, ■ ■ ■ , x^, t) for 
finding the set of eigenvalues near the values x\, - ■■ ,xn at time t, obeys the SFP 
equation 

The average density of eigenvalues, p(x), may be obtained from P by integrating 
over N - I variables. Specifically: 



r N N 

p(x,t)= I P[ dxu P(x\, ■ ■ ■ , xn, t) ^ 8(x 

*J i — 1 f_i 



N 

A 7 ). (99) 

k=l l=i 



with normalization J dxp(x) = N. Similarly we define the "two-particle" density 
p(x,y) = (Z/=i "Lj#i?>( x ~ x i)5(y - x j))> with J dxdy p(x,y) = N(N - 1). These 
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various densities obey an infinite hierarchy of equations obtained form Eq. ( [98] ) 
for P. Thus, the equation relating the one and two particle densities reads 

dp(x, t) 1 d 2 p(x, t) dC, p(x, y, t) 



f 



dy^^, (100) 



dt 2 dx 2 dx J x - y 

where j denotes the principal value of the integral. 

In the large N limit, this equation becomes a closed equation for the one parti- 
cle density To show that, we set 

p(x,y) = p(x)p(y) + p con (x,y), (101) 

where p C on(x,y) is the connected part of the two-point density. Then we change 
the normalization of the single particle density, defining 

p(x)=Np(x), (102) 

and similarly p(x,y) = N(N - l)p(x,y). At the same time, we rescale the time so 
that r = Nt |251. One then obtains 

n(x,y) 



dp{x) d £ p(y) 1 d 2 p(x) £ ^ Pc 
— + -p(x)j-^— = - — + 



x-y 2N dx 2 J x - y 

(103) 



In the large Af limit, the right hand side vanishes, leaving as announced a closed 
equation for pipe, r). Taking the Hilbert transform of the above equation, and fol- 
lowing the procedure from the previous sections, we immediately recognize the 
complex inviscid Burgers equation for the resolvent 



^'ii^l-f*-' <104) 

\N z-H(t)I J z-y 



which reads explicitly 



d T G(z, t) + G{z, r) d z G{z, r) = 0. (105) 

Let us comment the differences between this equation and the standard diffusion: 
First, the usual Laplace term (1/2) A has vanished, since in the matricial case this 
term is dwarfed by l/N factor. An additional non-linear term has however ap- 
peared, due to the interaction of diffusing eigenvalues, a term which by definition 
is absent in a one-dimensional random walk. This is how the inviscid (complex) 
Burgers equation [8l |25j appears in the matrix- valued diffusion process. The re- 
sulting nonlinearity can trigger shock waves, as we will see below. 
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Fig. 7. The Wigner semi-circle. 



Repeating the method of (complex) characteristics described in Sect. [3} with 
the characteristics determined by the implicit equation 

z = t + tG ($, Gq(z) = G(z, t = 0) = 1/z. (106) 

and assuming the solution £(z, t) to be known, the Burgers equation can be solved 
parametrically as G(z, r) = Gq(£(z, t)) = Go(z - tG(z, t)). The solution of this 
equation that is analytic in the lower half plane is 

G(z, t) - — (z - Vz 2 - 4r) , (107) 
2r 

whose imaginary part yields the familiar Wigner's semicircle for the average den- 
sity of eigenvalues. This is perhaps the fastest, and quite intuitive, derivation of 
this seminal result. 

In the fluid dynamical picture suggested by the Burgers equation, the edge of 
the spectrum corresponds to a singularity that is associated with the precursor of 
a shock wave, sometimes referred to as a "pre-shock" (23J. As discussed earlier, 
this singularity occurs when 

dz/d£ = = l+TG' (£ c ), (108) 
defining &(t). Since = &(t) = ±^, and 

z c = f c + TG (f c ) = ±2V?. (109) 

That is, the singularity occurs precisely at the edge of the spectrum, traveling with 
time t. Furthermore, the resulting singularity is of the square root type. To see 
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that, one expands the characteristic equation around the singular point. One gets 
Z - Zc = J(f - €efG^ c ) = J(f - Zc) 2 . (110) 

It follows that, in the vicinity of the positive edge of the spectrum z - z c = 2 y/r, 

^-^-±T 1/4 V^- (HI) 

Thus, as z moves towards z c and is bigger than z c , £ moves to £ c on the real axis. 
When z becomes smaller than z c , £ moves away from £ c along the imaginary axis. 
The imaginary part therefore exists for z < z c and yields a spectral density p(z) ~ 



s/z c - z, in agreement with ( 107 1. This square root behavior of the spectral density 
implies that in the vicinity of the edge of the spectrum, the number of eigenvalues 
in an interval of width s scales as Ns 3 ^ 2 = (N 2 ^s) 3 ^ 2 , implying that the interlevel 
spacing goes as N~ 2 ^. 

This scaling of the preshock wave resembles the Airy scaling in random matrix 
models. Indeed, we can provide a rigorous argument why this is the case |[24l . 
Note that in the case of an additive diffusion we can simply find the solution of 



Eq. (98 1 



l — r 

P(xi,---,x N ,f) = C\ UXi- xjfe' 1 '^, (112) 
i<j 

with C a (time-dependent) normalization constant. Indeed, in the random walk 
described above, the probability distribution retains its form at all instants of time. 
This means that we can repeat the standard stationary solution in terms of time 
dependent Hermite polynomials, which now remain orthogonal with respect to the 
time-dependent measure exp (-./Vx 2 /(2t)). Explicitly, the monic, time-dependent 
orthogonal Hermite polynomials read 

n k {x, r) = <=^ ^ J dq q k e-^\ (113) 



and satisfy 



r 



-— 2 

dxe I* n n {x,T)n m {x,T) = 6 nm c n , (114) 



with c 2 = n\ ^/2j:T n+[/2 , where we have used conventions from [26]. Note that the 
monic character of the /r„'s is not affected by the time dependence. 

By using the integral representation ( |1 13| >, it is easy to show that the n n (x, r)'s 
satisfy the following equation 

d T n n (x, t) = -v s d 2 x n n (x, t), (115) 
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with v s = 1 /2N, for any finite N. This is a diffusion equation with, however, a 
negative diffusion constant, which prevents an immediate physical interpretation. 
One can however understand intuitively this negative sign: a positive diffusion 
(viscosity) would smoothen the shock wave. In order to obtain the wildly oscillat- 
ing pattern of the preshock corresponding to the universal spectral fluctuations in 
the random matrix theory, we need an opposite mechanism. Note also that n n is 
an analytic function of x, and n„(-iy, r), with y real, satisfies a diffusion equation 
with a positive constant. 

As a last step, to see the emerging viscid Burgers structure, we perform a so- 
called inverse Cole-Hopf transform, i.e., we define the new function 

k 

Mz,t) = 2v s dMn k {z,r) = i V l — (116) 

N j-f Z~ Xi(T) 

with Imz + 0. The resulting equation for ft is the viscid Burgers [ 15 ] equation 

drfkiZ, T) + MZ, T)d-J k { Z , T) = -V s d 2 J k {z, T). (117) 



The equation ( 1 17 ) is satisfied by all the functions f k , for any k. We shall focus 
now on the function fy associated to 7In(z, t), due to the known fact that 7In(z, t) 
is equal to the average characteristic polynomial E9l . i.e 

<detfe-ff(T))>=*tf(z,T). (118) 

Note that in the large N limit, d z ln(det(z - H{t))) » <9 z (ln det(z - H{t))) = N G(z). 
Thus /n(z, t) coincides with the average resolvent G(z, r) in the large ,/V limit. In 



fact the structure of fy, as clear from Eq. ( 1 16 ), is very close to that of the resol- 
vent, with its poles given by the zeros of the characteristic polynomial. Eq. (117) 
for fy(z, t) is exact. The initial condition, fy(z, r = 0) = l/z, does not depend 
on N, so that all the finite ,/V corrections are taken into account by the viscous 
term. This observation allows us to recover celebrated Airy universality in the 
random matrix models, this time solely from the perspective of the theory of tur- 
bulence |[27l . Let us recall that in the vicinity of the edge of the spectrum, and in 
the inviscid limit, 

Mz,t)^±-^ + -^^Jz-=z~ c . (119) 



We set 



x = z c (j) + v 2 J 3 s, f N (x, t) = Zc(j) + v 1 J 3 xn(s, t), (120) 



with z c = d T z c = ± 1 / V^- The particular scaling of the coordinate is motivated by 
the fact that near the square root singularity the spacing between the eigenvalues 
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scales as N 2 ^ 3 . A simple calculation then yields the following equation for x(s, t) 
in the vicinity of z c (r) = 2 y/r: 



d T z c + v s ' ^-+Xw~ 

OT OS 



which, ignoring the small term of order v l J 3 , we can write as 



ds 2 ' 



d 

Is 



1 



dx 



' 2t 3/2 + 2 X + ds 



0. 



(121) 



(122) 



Note that the expression in the square brackets represents the Riccati equation, so 
the particular explicit solution corresponding to characteristic polynomial is 



X(s, t) = 2 



Ai'(a 1/3 (j)) 
Ai{a l ' 3 {s)) ' 



(123) 



where Ai denotes the Airy function, and a = 1/(4t) 3 ^ 2 . 

For completeness we note, that the Cauchy transforms of the monic orthogonal 
polynomials 



Pk(z 



t) = — r 

2m J_ 



7Tft(jC, r)e 



-Nx 2 /2t 



x-z 



(124) 



which are related to average inverse spectral determinant, generate similar univer- 
sal preshock phenomena in complex Burgers equations. Their detailed behavior in 
the vicinity of the preshock is related to the two remaining solutions of the Airy 
equation (Airy of the second kind, [28]). 



5.2. Multiplicative random walk of large matrices 

The case of multiplicative (geometric) one dimensional random walk can be 
easily reduced to the additive one, since e x+y = e x e y , so taking the logarithm of the 
multiplicative random walk reduces this case to the additive random walk of the 
logarithms of multiplicative increments. Note that this simple prescription does 
not work in the case of matrices. First, the product of two hermitian matrices is no 
longer hermitian, second, matrices do not commute, so e x+Y + e x e Y . This means 
that the matrix-valued geometric random walk may exhibit new and interesting 
phenomena. To avoid the nonhermiticity of the product of hermitian matrices, we 
stick to unitary matrices - their product is unitary and a simple realization is 

U - exp{i yfdtH] (125) 

where H is hermitian. Note that this procedure is equivalent to Janik-Wieczorek 
model, described in the first part of these notes. Following Dyson, we recover that 
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the only difference corresponding to the multiplicative case is the modification 
of the electric force. Since the eigenvalues of the unitary matrices are forced to 
stay on the unit circle, the electric force has to respect periodicity, i.e. has to 
gain contributions from all the distances between the interacting charges modulo 
2nm, where m is the integer. This infinite resummation modifies the electric force, 
yielding a SFP equation for Circular Unitary Ensemble [21] and corresponding to 
the process where 

(6y i ) = E(y i )At, ((8yd 2 ) = At, (126) 

where E(yi) = ^ T*i?j c °t(yi ~ yj)/2, parameterizing the eigenvalues on the unit 
circle as z = exp((y)- This particular form of the electric field was already met 
in the first part of these lectures. Let us make two remarks. First, we note that 
in the present case, in order to reach the equilibrium solution, corresponding to 
a uniform distribution of eigenvalues along the unit circle, one does not need the 
stabilizing external force - due to the compactness of the support of the eigenvalues 
the mutual repulsion between them is sufficient to reach equilibrium in the infinite 
time limit. Second, we may expect a new universal phenomenon - due to the 
compactness of the support of the spectrum, left and right Airy preshocks at the 
ends of evolving packet of eigenvalues have to meet at some critical time, creating 
universal generalized Airy function exhibiting novel scaling with large N. 

In principle, we can follow the route sketched above for the additive case, 
i.e. find the time-dependent solution of the pertinent SFP equation and construct 
orthogonal polynomials. This time however the mathematics is more involved - 
the time dependent solution does not factorize nicely, and depends on the ratio of 
Vandermonde determinants built from Jacobi theta functions (solutions of diffusion 
equation with periodic boundary conditions), and the corresponding orthogonal 
polynomials belong to Schur class. In order to avoid these complications and to 
make connection with the lectures of H. Neuberger |[30l at this School, we make 
the following observation [3TJ . Let us define a generic observable 




(127) 



where the sum is over the representations, xr denotes the character of the repre- 
sentation R and djt its dimension. Calculating the expectation value of the above 
object with respect to the measure i30l 

P(U) = ^ d RXR (U)e~ TC ^ R \ (128) 

R 

where C2(R) is a Casimir operator for the representation R, we get 

(e- aT 0(z, U)) = e-" T ^ a R (z) e~ TC ^ R \ (129) 

R 
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If we act now on this expression with the heat-kernel operator (zd z ) 2 + k8 t , we 
notice that the heat equation does not mix the coefficients cir corresponding to 
different representations. This allows us to write down a simple solution 



This expression has to be well defined on the unit circle, which implies that the 
argument of the square root has to be the square of an integer. This imposes 
very strong restriction on Young diagrams forming the representation R. For an- 
tisymmetric representations considered in Ref. IfTTll . the single valued- ness of the 
solution corresponds to a = -N/8, provided the evolution time is identified as 
t = t(l + l/N). Indeed, the combination of Casimir and a reads in this case 



Jn \1 ~ ^) ' f° r an y ^ runn i n 8 between and N. This gives the whole family of 
(N + 1) heat equations, one for each k, in analogy to the family of monic Her- 
mite polynomials fulfilling the heat equation for the case of additive random walk. 
Since the heat equation is linear, any combination of monomials in z fulfills the 
same equation. In particular this is the case of the characteristic polynomial con- 
sidered by Neuberger IfTTll . As before, the inverse Cole-Hopf transformation yields 
the viscid Burgers equation. Note that this construction explains naturally the ap- 
pearance of the somewhat mysterious extra factor e~ Nr ^ in the corresponding vis- 
cid Burgers equations for the characteristic polynomial. For completeness we note 
that a similar viscid Burgers structure appears for the inverse characteristic poly- 
nomial |[T8l . Now, the relevant representations are symmetric, which corresponds 
to "time" t = ?(1 — jj), and unrestricted k running from zero to infinity, and 
the argument of the square root is proportional to (j + k) 2 , provided that again 
a = -N/8. This gives the infinite hierarchy of heat equations for any k. Since 
1 / det(z - U) can be expanded into infinite series in 1 /z, again the pertinent inverse 
Cole-Hopf transformation reproduces viscid complex Burgers equation, as noted 
by Neuberger. 

Alike in the additive case, we can now analyze the critical behavior at the 
closure of the spectral gap, using the tools of the theory of turbulence for the coa- 
lescence of two Airy-like universal preshocks merging at z = - 1 at some critical 
time. However, instead of going into the beautiful, but involved mathematical con- 
struction corresponding to the appearance of the universal Pearcey function at the 
closure of the gap, we propose to exploit a somewhat unexpected analogy between 
the spectral properties of the gap in large N c Yang-Mills theory and geometric 
optics. 

6. Diffraction catastrophes and large N universalities in YM theories 

The large N limit is very often considered as a "classical" one, since it switches 
off fluctuations. In the spirit of this analogy we consider classical (geometric) op- 
tics, where the wavelength vanishes, A = 0, and rays of light are straight lines. 



.±TjK(C 2 (R)+a) 



(130) 
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Type 








Fold 


P=i 




iff = 2nAi(g) = J_ m dt exp «'(f 3 /3 + £0 


Cusp 




°~X=\ O-y = f 


</r = />(£ 77) = jT A exp i(t 4 /t + f?/2 + rjt) 



Table 1. Classification of two lowest stable singularities 



These rays of light can condense on some surfaces, yielding high intensity (actu- 
ally, in the A = limit infinite intensity) hypersurfaces, which are called caustics, 
from a Greek word denoting "burning". If we relax the restrictions of geometric 
optics, i.e., if we allow for some very small A + 0, wave patterns of light appear. 
Rays start to interfere, forming wave packets, and intensity is no longer infinite. 
We may ask now, how the limiting procedure from wave optics to geometric optics 
takes place. In other words, we ask how the wave packet scales with A — > 0. The 
answer to this intriguing question has been provided by Berry and Howls (32), 
and their classification of "color diffraction catastrophes" is a particular, physical 
realization of the classification of singularities in the so-called catastrophe the- 
ory. There are only seven stable caustics, the two lowest, relevant for our analy- 
sis, corresponding to so-called fold and cusp singularities, controlled by one and 
two parameters, respectively . The scaling of the wave packet in these two cases 
reads (351 



(131) 



where r denotes a set of "control parameters" and the universal critical exponents 
<Ti and /3 are known respectively as Berry and Arnold indices. The Table 1 . summa- 
rizes the universal properties of caustics. Immediately we see an analogy between 
the Airy (Ai(g)) and Pearcey (P(£, rf)) universal functions appearing at the fold and 
cusp ( merging of two folds) and which describe the properties of the character- 
istic polynomials for critical spectra of Wilson loops in Yang-Mills theory, also 
for D > 2. Complex characteristics (straight lines) play the role of rays of line. 
Singularities (shock waves) correspond to caustics, finite viscosity (1/2N) scaling 
has the same critical indices as finite A scaling for the intensity of the wave pack- 



ets, see Eqs.( 82 )-( 84 ). The exact correspondence is summarized in the form of the 
Table 2. 

We find it rather intriguing that the beautiful scaling of interference fringes 
in diffractive optics could belong to the same universality class as the finite N 
critical scaling of the spectral density of the Wilson operator in non-Abelian gauge 
theories in nontrivial dimensions. 
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GEOMETRIC OPTICS 


jV - oo YANG-MILLS 


wavelength A = 


viscosity v s - - 


rays of light 
caustics 


rays of characteristics 
singularities of Wilson loop spectra 


WAVE OPTICS 


FINITE N YANG-MILLS 




< det(z - W(Q) > 


Fold scaling cr = 2/3 


N 2/i scaling at the edge 


Cusp scaling cr\ = 1/2 <T2 - 3/4 


N 112 and N i/4 scaling at critical size 



Table 2. Morphology of singularities - analogies. 




Fig. 8. A familiar example of caustics Il33l . Note two fold-line caustics merging into the 
cusp. 

7. Outlook 

In these lectures, we have shown how the study of the complex Burgers equa- 
tion sheds a new light on universal properties of the the large N c transition in multi- 
color Yang-Mills, that was first identified by Durhuus and Olesen many years ago. 
These lectures should be viewed as complementary to those by Neuberger and by 
Narayanan at this School, and they offer a view of a similar subject from differ- 
ent angles and perspectives. Our observations allow us to link together domains 
of theoretical physics and mathematics that are rather unrelated at first sight: The 
spectral flow of eigenvalues of the Wilson operators has features reminiscent of 
classical turbulence, the universal behavior is locked by multiplicative unitary dif- 
fusion considered by Dyson already in 1962, critical exponents belong to a classi- 
fication of stable singularities (catastrophe theory) and the phenomenon of scaling 
with finite Af at critical size of the Wilson loop has exact counterpart in diffractive 
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Fig. 9. Color diffraction catastrophe ||34ll . Picture shows modulus of Pearcey function, 
anisotropic character of interference fringes reflects different scaling exponents cr\ and 

optics! The fact the whole dynamics of complicated non-perturbative QCD can 
be reduced in some spectral regime to a matrix model is not new - a notable case 
is the universal scaling of the spectral density of the Euclidean Dirac operator for 
sufficiently small eigenvalues, where the spectrum belongs to the broad universal- 
ity class of the corresponding chiral models 11361 . What we find remarkable is that, 
in the case of the large N c transition, the analogous universal matrix model seems 
to be represented by the simplest realization of the multiplicative matrix random 
walk. In general, one may expect that in a very narrow spectral window around 
A = - 1 a universal oscillatory behavior appears as a preparation for the formation 
of the spectral shock-wave, in qualitative analogy to similar spectral oscillations 
of quark condensate before spontaneous breakdown of the chiral symmetry based 
on Banks-Casher relation ll37ll . 

Certainly, further generalizations and analogies are possible. We mention here 
the extension to supersymmetric models, intriguing role of the fermions [38 ] and 
the analogies with shock phenomena in mezoscopic systems (universal conduc- 
tance fluctuations) [39 ] or growth processes of the Kardar-Parisi-Zhang universal- 
ity class and statistical properties of the equilibrium shapes of crystals [40]. We 
hope that these lectures will trigger the need of better understanding of all these 
relations and analogies. 
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